Statistics of the Number of Zero Crossings : from Random Polynomials to Diffusion 

Equation. 
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We consider a class of real random polynomials, indexed by an integer d, of large degree n 
and focus on the number of real roots of such random polynomials. The probability that such 
polynomials have no real root in the interval [0, 1] decays as a power law n~ e ^ where 6(d) > is 
the exponent associated to the decay of the persistence probability for the diffusion equation with 
random initial conditions in space dimension d. For n even, the probability that such polynomials 
have no root on the full real axis decays as n _2 ^ £, ' d ' +e ' 2 ^ . For d = 1, this connection allows for a 
physical realization of real random polynomials. We further show that the probability that such 
polynomials have exactly k real roots in [0,1] has an unusual scaling form given by ti _| pw logn ) 
where (p(x) is a universal large deviation function. 

PACS numbers: 02.50.-r, 05.40.-a,05.70.Ln, 82.40.Bj 

Persistence properties and related first passage prob- 
lems have been the subject of intense activities, both 
theoretically [l[ and experimentally 0, Q these last few 
years. The persistence probability p(t) for a time depen- 
dent stochastic process of zero mean, is defined as the 
probability that it has not changed sign up to time t. In 
many physical situations, p(t) was found to decay at large 
time as a power law p(t) cx t~ e . Surprisingly, the persis- 
tence exponent 9 was found to be highly non trivial even 
for simple systems. One example is the diffusion equation 
in space dimension d, dt<t>(x,t) = S7 2 cj)(x,t) with "white 
noise" initial conditions (4>{x, 0)(p(x' , 0)) = S d (x — x'). 
The solution is characterized by a single growing length 
scale £(t) cx t x l z , with z = 2. For a system of linear size 
L, the persistence po(t, L) is simply the probability that 
(Mx,t) does not change sign up to time t. It was found 
[J] that for t ^> 1, po(t, L) has a finite size scaling form 



physical realization of the complex roots of Weyl complex 
polynomials in a system of a rotating quasi-ideal atomic 
Bose gas. Here we focus instead on the real roots of a 
class of real random polynomials indexed by an integer d 



f n (x) = a + ^ a i 



(2) 



p Q {t,L)<xL- ze ^h(L z /t) 



(1) 



with h(u) ~ c st , a constant independent of L and t, for 
u -C 1 and h(u) cx u e ^ for u 1, where 0(d) is a non 
trivial exponent, e.g. 0(1) = 0.1207,0(2) = 0.1875. The 
scaling form of po(t, L) indicates that pa(t, L) cx t~ e ^ for 
large t in an infinite system. Alternatively, 0(d) can be 
extracted by measuring po(t, L) cx L~ ze ( d ) for very large 
time t ^> L z . Remarkably, the persistence for d = 1 
was observed in experiments on magnetization of spin 
polarized Xe gas and the exponent 9 eKp (l) — 0.12 was 
measured Q , in good agreement with analytical approx- 
imations and numerical simulations 0, [1| . 

Another apparently unrelated problem concerns the 
roots of random polynomials (i.e. polynomials with ran- 
dom coefficients), which have attracted renewed interest 
over the last few years 0, 0] in the context of proba- 
bility and number theory. A recent work 0] proposed a 



Here a^s are real Gaussian random variables of zero mean 
and with correlations (aidj) — Sij. We will see below 
that, for d = 2, where f n (x) reduce to the famous Kac 
polynomials [9( , the statistics of the real roots of /„ (x) is 
identical in the 4 subintervals ] — oo, — 1[, [— 1, 0], [0, 1] and 
]1, +oo[. However, for d ^ 2, the statistics of real roots of 
f n (x) depend on d in the two inner intervals [—1,0] and 
[0, 1], while it is identical to the case d = 2 in the two 
outer ones. In this letter we will focus primarily on the 
interval [0, 1] and ask : what is the probability Po(l,n) 
that f n {x) has no real root in [0, 1] ? Recently, it was 
found, for d — 2, that Po(l, n) oc n~^ 2 ) for large n where 
the exponent £(2) ~ 0.19(1) was computed numerically 
[lfl ]. In addition, for the special case d = 2, the authors 
of Ref. [l(| showed that Po(lj n ) is related to the proba- 
bility of no zero crossing of a Gaussian stationary process 
(GSP) with correlator sech(|T|/2). 

The purpose of this Letter is to provide a link between 
the persistence of the diffusion equation and the proba- 
bility Pq(1, n) that f n (x) has no real root in [0, 1]. For ar- 
bitrary dimension d, we show that Po(l, n) cx n~ ( ' l - d > with 
((d) = 0(d). Given that 0(1) was measured experimen- 
tally [![, this demonstrates an experimental realization 
of real random polynomials. The connection between 
these two problems in arbitrary d is achieved by showing 
that both problems can be mapped to the same GSP. In 
addition, we compute the probability that a "smooth" 
GSP, such as the one that appears in the context of dif- 



2 



fusion equation, crosses zero exactly k times up to time 
T . Translated into the language of random polynomi- 
als, our analysis shows that the probability Pfc(l, n) that 
f n (x) has exactly k real roots in the interval [0,1] has 
a rather unusual scaling form (for large k, large n, but 
keeping the ratio k / log n fixed) 



P fe (l, n) oc n 



log 71 



(3) 



where <p(x) is a large deviation function, with ^3(0) = 
C(d). Besides, our numerical analysis suggests that <p(x) 
is universal in the sense that it is independent of the 
distribution of a,i provided (a 2 ) is finite. 

To study the persistence probability po(t, L) of the dif- 
fusion equation, it is customary to study the normalized 
process X(t) — c/)(x 1 t)/((j)(x,t) 2 ) 1 / 2 [I]. Its autocorrela- 
tion function a(t,t') = (X(t)X(t')) is given, in the limit 
t,t' < L 2 by a(t,t') = [Att'/(t + t') 2 ] d / i . In terms of loga- 
rithmic time variable T = logt, X(T) is a GSP with cor- 
relator a (T, T) = a(T-T') = [sech(|T-T'|/2)] d / 2 , which 
decays exponentially for large \T — T'\. Thus the persis- 
tence probability p n (t,L), for t <C L 2 , reduces to the 
computation of the probability Vo(T) of no zero cross- 
ing of X(T) in the interval [0,T]. It is well known [Tl | 
that if a(T) < l/T at large T then V (T) ~ exp[-6»T] 
decays exponentially for large T where the decay con- 
stant 9 depends on the full stationary correlator a(T). 
Reverting back to the original time t = e T , one finds 
po(t,L) ~ t~ e{ - d \ for t < L 2 . In the opposite limit 
t 3> L 2 , one has po(t,L) — > Aj,, a constant which de- 
pends on L. These two limiting behaviors of p$(t,L) 
can be combined into a single finite size scaling form in 
Eq. (P) where 9(d) is the decay constant associated with 
the no zero crossing probability of the GSP with correla- 
tor a(T) = [sech(|T|/2)] d / 2 0|. 

The mapping of real random polynomials f n (%) in 
Eq. j2} to a GSP is more subtle. We first observe 
that for large n the real roots of f n (x) are concentrated 
around x = ±1. To show this, we have generalized the 
Kac's method [9] to compute the mean density of real 
roots p n (x) of f n {x) ([2]). We do not give the details 
of the computation and simply quote the results here. 
We find that p„(±l) diverges in the large n limit as 
p n (±l) ~ 2n-\/ dj (d + 4)/ (ir(d+ 2)) and away from these 
singularities Poo(x 7^ ±1) is given, for \x\ < 1, by 



Poo(x) 



(U^ d/2 (x 2 )(l + U^jx 2 )) - U 2 _ d/2 (x 2 ))2 
ir\x\(l + Lii_ d / 2 (x 2 )) 

(4) 

where Li„(z) = y^ —1 z l ji n is the polylogarithm function, 
yielding back p^ix 7^ ±1) = (tt\1 — s 2 |) _1 for d = 2. In 
particular, one has /9oo(0) = 1/tt for all d, and /foo^) ~ 

(d/2)2(2 7 r((l-a:)))- 1 fbrx-^l". For \x\ > 1, Poc {x) = 
1/[tt(x 2 - 1)] for all d [13]. In Fig. Q]a), we plot poo(a;) 
for d = 1 where the divergence at x = ±1 indicates that 
the real roots concentrate around x = ±1 for large n. 
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FIG. 1: a) : poo(x) given analytically in Eq. Q and below 
as a function of x for d = 1. The divergences for x = ±1 
indicate that the real roots concentrate around a; = ±1 for 
large n. b) : Plot of Po(l,n) as a function of n for different 
values of d — 1, 2, 3, 4, 10. The lines are guide to the eyes. The 
measured exponents are given by 0.12, 0.18, 0.23, 0.27, 0.46 for 
d = 1,2,3,4,10 respectively, in full agreement with the nu- 
merical values for 6(d) obtained in Ref. [J, [g. 



The random polynomial f n (x) being a Gaussian pro- 
cess, it is completely determined by its two-point corre- 
lator C n (x,x') — (f n (x)f n (x')). Using the fact that Oj's 
are uncorrelated, it is easy to see that 



C n (x, x ) 



E 



(d-2)/2 



(xx'y. 



(5) 



It is useful to introduce the normalized correlator 

1 1 

Cn(x, X ) = C n (x,x')/C n (x,x)2C n (x' ,x')z . From the 
analysis of C n (x, x 1 ) in the large n limit, one can show 
that, asymptotically, f n (x) takes independent values in 
the 4 subintervals ] — 00,— 1[, [—1,0], [0,1] and ]l,+oo[. 
For d — 2, C n (x, x') has a special symmetry, namely 
C n (x,x') — C n (—x,—x') = C n (l/x,l/x'), meaning that 
the normalized Gaussian processes in these 4 subinter- 
vals are independent and isomorphic. Denoting P (l,n) 
as the probability of no zero crossing in [0,1], it follows 
that the probability go ( n ) that there are no real roots for 
d = 2 is precisely equal to [i- u (T"')] 4 for large n. For 
d 7^ 2, one can show that the behavior of C n (x,x') de- 
pends on d in the two inner intervals (see below), while 
it behaves like for d = 2 in the two outer ones [16( . 

We now focus on the interval [0, 1] and to make a pre- 
cise connnection with the persistence probability po (t, L) 
defined in the context of the diffusion equation we define 
P (x,n), for x < 1, as the probability that f n (x) has no 
real root in the interval [0, x]. We next reparametrize 
the polynomial with a change of variable, x = 1 — 1/t. 
One finds that the relevant scaling limit of C n (t, t') is ob- 
tained for t,t',n — ► 00 keeping t = t/n and r = t'/n 
fixed (see also Ref. [13] )■ I n ih&t scaling limit one finds 
that C n {t,t') — > C(t,i') with the asymptotic behaviors 



C(t,t') 




1 1' < 1 

t, V > 1 



(6) 



For t, t' <C 1, this correlator is exactly the same as the 
one found for diffusion, C(t,t') = a(t,t') in this regime. 
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FIG. 2: Plot of n i2) P {x,n), with 0(2) = 0.1875 as a function 
of (1 - z)n for different degrees n = 64, 128, 256, 512, 1024. 
Inset : Po(x,n) as a function of 1 — a; for different degrees n. 



Since a Gaussian process is completely characterized by 
its two point correlator, we conclude that the diffusion 
process and the random polynomial are essentially the 
same Gaussian process and hence have the same zero 
crossing properties. In the opposite limit, t,r 1 the 
fact that C(t, r) — > 1 suggests that Po(x,n) goes to a 
constant when x — > 1. Therefore, in complete analogy 
with Eq. |T]) we propose the scaling form for random 
polynomials 



P (x,n) oc n- e{d) h(n(l - 
" s * for u <C 1 and 



(7) 



with /i(u) ~ c st for u <g; 1 and h(u) ~ for u ^> 1, 

where 8(d) is the persitence exponent associated to the 
diffusion equation in dimension d. Note that n here plays 
the role of L z in diffusion problem while the variable 1 — x 
is the analogue of the inverse time l/t. This scaling form 
([7]) , which we verify numerically, thus establishes a direct 
link between the real roots of real random polynomials 
and the diffusion equation with random initial conditions. 
It also follows from our earlier discussion on the 4 subin- 
tervals that the probability that f n (x) has no real root 
decays like n" 6 ^) with b(d) = 2(8(d) +0(2)) where 26(d) 
is contributed by the 2 inner intervals and 29(2) by the 
two outer intervals. 

We have verified this scaling form ([7]) numerically by 
computing the number of real roots in the interval [0, x] 
of random polynomials such as in Eq. @ for differ- 
ent degrees n. In each case, the probability distribution 
Pq(x, n) is obtained by averaging over 10 realizations of 
the random variables a^'s, drawn independently from a 
Gaussian distribution of unit variance. The plot shown 
on Fig. [2] for d = 2, shows a good agreement with the 
scaling in Eq. ((Jj with 6(2) — 0.1875, in agreement with 
the numerical value reported in Ref. [Icj . We have also 
calculated the exponent 0(d) in Eq. (jTJ) , for various di- 
mensions, by measuring Pq(1, n) for different values of d. 



The results are shown in Fig. [T] b) where, according to 
our prediction in Eq. ([7]) Pq(1,tl) ~ n~ e< * d ) for large n. 
The value of 8(d) obtained this way is fully compatible 
with previous numerical estimates H H of 8(d) for the 
diffusion equation. 

We now generalize our analysis and consider the prob- 
ability Pk(t, L) that the diffusing field 4>(x, t) crosses zero 
exactly k times up to time t (similarly, one considers the 
probability Pk(x,n), x < 1, that such polynomials (J2J> 
have exactly k real roots [1(3] in the interval [0, x], see be- 
low) . Let us first consider the regime 1<(<I Z . In this 
regime, pk(t,L) is given by the probability Vk(T) that 
X(T) crosses zero exactly k times where X(T) is a GSP 
with correlations a(\T-T'\) = [sech(|T - T'\)] d / 2 , where 
T = logt. Since, a(T) = 1 - ^T 2 + o(T 2 ) for small T, 

our GSP is a smooth process with a finite density of zero 

l 

crossings given by the Rice's formula fj, = (— a"(0))2 /it 
(l3t - We propose the following scaling form for large T 
and large k 



\ogV k (T) = -Tip 



k 

7? 



(8) 



To understand the origin of this scaling form, let us con- 
sider the generating function V(p,T) — '}2 k K ' =0 p k 'P(k,T) 
[l3| . It turns out that, V(p,T) - exp(-%)T), where 
for a smooth GSP 8(p) depends continuously on p : this 
was shown exactly for the random acceleration process 
and approximately using the independent interval ap- 
proximation for arbitrary smooth GSP - and checked 
numerically for the diffusion equation with random ini- 
tial conditions 



14|. 



If the scaling in Eq. (|8|) holds, 
one gets by steepest descent method valid for large T, 
9(p) = Min a;> o[/ia: logp — (fi(x)]. Inverting the Legendre 
transform we get 



ip(x) = Max <p< 2 [AM;logp + 8(p)] 



(9) 



Notice that although 9(p) is a priori defined on the inter- 
val [0, 1], the computation of ip(x) involves an analytical 
continuation of 9(p) on [0, 2]. Going back to real time t, 
Eq. © then yields a rather unusual scaling form valid 
in the limit 1 < t < L 2 



Pk(t,L)~t 



' ( A» log t ) 



(10) 



For k close to [i\ogt, one expects pk(t,L) to behave lo- 
cally as a Gaussian and ip(x) is thus quadratic around 
x = 1. Away from the minimum, we have not been 
able to obtain ip(x) analytically. We have thus tested 
the scaling form (fT0|) numerically. We used a space time 
discretized diffusion equation 



4>i (t + 1) = fa (t) + a [<t>j (*) ~ <t>i (*)] 

3 



(11) 



where j runs over the nearest neighbours of i on a d- 
dimensional square lattice. A stability analysis shows 
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that the solution is unstable for a > a c = l/2d and we 
chose a — a c /2 that provided the quickest onset of the 
asymptotic behavior. The initial values </>j(0)'s were cho- 
sen independently from a Gaussian distribution of zero 
mean and unit variance. Fig. [3J shows the results for 
d = 2, i.e. the case of Kac's polynomials ©. We com- 
puted Pk(t, L) for different times f for a system of linear 
size L = 256 and averaging over 100 different realizations 
of the random initial condition. In the inset of Fig. [3j 
we plot (— logpk(t, L))/ logf as a function k where the 
different curves correspond to different times f . In Fig. [3J 




1 2 3 4 5 6 7 
2 jtk/log(t) 

FIG. 3: (— logpfc(t, L))/ logt for the diffusion equation with 
random initial conditions as a function of fc//ilogi with 
H = l/2ir for different times f = 256, 512, 1024, 2048. Inset : 
(— \ogpk(t, L))/ logf as a function of k for different times. 

these curves for different time f fall on the same master 
curve when plotted as a function of the rescaled variable 
k/plogt with fi = 1/27T in this case, confirming the va- 
lidity of the scaling form in Eq. (flC)]) . We also checked 
that a different choice of the distribution of <fii(0), such 
as 4>i(0) = ±1 or rectangular, gave, within the error bars, 
the same function if(x) thus indicating the universality 
of this large deviation function. 

In the opposite limit f > L z , one simply replaces f in 
(|10P by L z . Translating into random polynomials, this 
regime corresponds to (1 — x) -C nT x since one just re- 
places f by 1/(1 — x) and L z by the degree n as dis- 
cussed before. Thus, in this regime, we arrive at the 
announced scaling form for qk(n) in Eq. ([3]). For the 
special case of Kac's polynomials (d = 2), this scaling 
form, in the neighbourhood of k = logn/27r, is consis- 
tent with the rigourous result [151 ] that in this neigh- 
bourhood qk(n) is a Gaussian with mean logn/27r and 
variance V n ~ — (1 — — ) logn in the large n limit. 

In fact, the scaling in Eq. (fTO]) holds more generally 
for any smooth GSP. To illustrate this, we considered the 
random acceleration process d 2 x(t) / dt 2 = 77(f) where 77(f) 



is a white noise, and for which p = y/3/(2n). For this 
particular smooth GSP 9(p) has been computed exactly 
[Ht, yielding §(p) = |(1 - | s in _1 (§)). By performing 
the Legendre transform ([9]) one obtains the asymptotic 
behaviors as 

(l + ^xlogx , x^0 

^)~|lgf(*-l) 2 > s-1 (12) 
[^2; log 2 , x — > 00 

which gives back the exact result (p(0) = 1/4 [l9j . 

To conclude, we have established in this Letter a con- 
nection between the persistence probability po(t, L) for 
the diffusion equation in dimension d and the probability 
Po{x, n) that generalized Kac's real random polynomials, 
indexed by d as in Eq. l|2|). have no real roots in the in- 
terval [0, a;], with x < 1. This connection is useful in pre- 
dicting new results for random polynomials, such as the 
unusual scaling form ([3J) for the probability of having k 
real roots in [0, 1]. Besides, we hope that this connection 
may also shed some light in calculating the exponents 
8(d) exactly which still remains a challenge. 

S.N.M thanks J. Unterberger for pointing out Ref. flfl ]. 
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